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Abstract. Using a hydrodynamic theory that incorporates a momentum decay mechanism, we calculate 
the drag coefficient of a circular liquid domain of finite viscosity moving in a two-dimensional membrane. 
We derive an analytical expression for the drag coefficient which covers the whole range of domain sizes. 
Several limiting expressions are discussed. The obtained drag coefficient decreases as the domain viscosity 
becomes smaller with respect to the outer membrane viscosity. This is because the flow induced in the 
domain acts to transport the fluid in the surrounding matrix more efficiently. 

PACS. 87.16.D- Membranes, bilayers and vesicles - 87. 16. dp Transport, including channels, pores, and 
lateral diffusion - 68.05.-n Liquid-liquid interfaces - 66.20.-d Viscosity of liquids; diffusive momentum 
transport 



1 Introduction 

Membrane components such as lipids and proteins are 
subject to Brownian motion, and the resulting diffusive 
process is an important mechanism for the transport of 
these materials. Consequently, it is generally believed that 
many biochemical functions in membranes are diffusion 
controlled processes. Up to now, the studies of diffusion 
in biomembranes have mainly concentrated either on pro- 
tein molecules 11121,31115] or, more recently, on lipid do- 
mains sometimes called as "lipid rafts" |6I7I8I9I1U] , Lipid 
domains appear due to a lateral phase separation between 

' a liquid-ordered phase and a liquid-disordered phase |llj , 
and have attracted much attention in relation to biological 

■ functionalities [T2"] . 

In the theoretical studies, a protein molecule is as- 
sumed to be a rigid disk moving in a two-dimensional (2D) 
thin fluid sheet under low Reynolds number conditions. 
Saffman and Delbriick (SD) investigated this problem by 
considering a membrane that is sandwiched in between a 
three-dimensional (3D) fluid medium [13114) . They showed 
that the drag coefficient has only a weak logarithmic de- 
pendence on the disk radius in the small size limit. The 
SD model was further analyzed by Hughes et al. for the 
whole range of protein sizes, and they also obtained the 
asymptotic expression in the large size limit |15] . The SD 
model was extended by Lcvine et al. for a viscoelastic film 
within an infinite viscous liquid, taking into account both 
in-plane and out-of-plane dynamics |16I17I18) . 

a E-mail: komura@tmu.ac.jp 



On the other hand, Evans and Sackmann looked at 
a slightly different situation [19], i.e., the diffusion of a 
protein molecule moving in a supported membrane (in- 
stead of a free membrane). The presence of the solid sub- 
strate is taken into account through a friction term in 
the Stokes equation. The equivalent hydrodynamic model 
was independently proposed by Izuyama who suggested 
that the momentum decay mechanism should generally 
exist for membranes surrounded by water |20) . Since a hy- 
drodynamic equation with a momentum decay term was 
originally proposed by Brinkman for flow in porous me- 
dia [2T] , we call these approaches as the 2D Brinkman- type 
model. One of the advantages of the Brinkman model is 
that the drag coefficient can be analytically obtained over 
the whole range of disk sizes. In the small size limit, the 
SD model and the Brinkman model give the same loga- 
rithmic dependence. In biological systems, the Brinkman 
model can be more relevant because the cell membranes 
are strongly anchored to the underlying cytoskeleton, or 
are tightly adhered to other cells |22) . 

Most of the previous work has concentrated on the 
problem of protein diffusion in a membrane for which the 
assumption of a rigid disk is reasonable. Liquid lipid do- 
mains, on the other hand, cannot be considered as rigid 
objects. A proper description should take into account the 
fluid nature and finite viscosity of the diffusing domain. 
It should be also noted that the viscosity of the lipid do- 
main is different from that of the matrix [7]. In this paper, 
we derive the drag coefficient of a circular liquid domain 
by taking into account its finite viscosity using the 2D 
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Brinkman-type hydrodynamic equation. We show that the 
drag coefficient decreases as the domain viscosity becomes 
smaller. An analogous problem for a spherical drop of fluid 
was solved by Rybczyhski and Hadamard about 1 00 years 
ago [23124) . Some authors considered the hydrodynamics 
of Langmuir monolayer domains at the air/water inter- 
face by taking into account the fluid nature of the do- 
main [25 26 27]. In these works, however, the domain and 
the matrix viscosities are assumed to be the same, and the 
drag coefficient was not obtained [28]. Although the rota- 
tional diffusion coefficient was also calculated for a rigid 
disk [13114115119] . such a motion is irrelevant for a fluid 
domain which has an internal flow. 

In the next section, we begin with a description of 
the 2D Brinkman model. In sect. [21 the solutions for the 
inside domain are connected to those for the outer mem- 
brane matrix through the appropriate boundary condi- 
tions. Then the drag coefficient is calculated for the entire 
region of the domain sizes. We also obtain several limit- 
ing expressions; some of which coincide with the previous 
results. Several discussions are provided in sect. |H 



2 Hydrodynamic model with momentum 
decay 

The usual method of obtaining the diffusion coefficient D 
from the drag coefficient £ is to use the Einstein relation 



D 



k B T 

C ' 



(1) 



where fce is the Boltzmann constant and T the tempera- 
ture. However, for a pure 2D system in the low Reynolds 
number regime, a linear relation between the velocity and 
drag force cannot be obtained [25J. This fact is due to 
the inability to simultaneously satisfy the boundary con- 
ditions both at the disk (cylinder) surface and at infinity 
within the Stokes approximation [30] . Such a problem is 
known as the Stokes paradox which essentially originates 
from the constraint of momentum conservation in a pure 
2D fluid. 

Although a lipid bilayer membrane can be treated as a 
2D viscous fluid, it is not an isolated system because the 
membrane is surrounded by a 3D fluid. Due to the cou- 
pling between the 2D membrane and the 3D fluid, the mo- 
mentum within the membrane can leak away to the outer 
fluid. Such an effect can be taken into account through 
a momentum decay term in the equation of motion. This 
also enables the 3D fluid to be discarded and we can work 
in a pure 2D system. A hydrodynamic equation which is 
consistent with the total momentum decay is 



<9v 



V)v 



ryV v — Vp — Av. 



(2) 



In the above, v and p are the 2D velocity and pressure, 
respectively, while p is the 2D membrane density, r\ the 
membrane 2D viscosity, and A the momentum decay pa- 
rameter. Notice that the usually reported membrane 3D 



viscosity is given by T]/h, where h is the membrane thick- 
ness. A more detailed derivation of eq. ([2]) in terms of 
the total momentum decay is given in [31[ . Since 3D ver- 
sion of eq. (0) was originally used by Brinkman for flows 
in porous media |21) . we call it as 2D Brinkman equa- 
tion. But it should be kept in mind that there is no real 
connection between porous media and membranes except 
from the momentum decay mechanism. For a typical flow 
in a membrane, we can adopt the Stokes approximation 
to eq. ([2]). Then we have 



7?V 2 v - Vp - Av = 0. 



(3) 



On the other hand, the incompressibility condition is given 

by 

V ■ v = 0. (4) 

Due to the presence of the momentum decay term, the 
Stokes paradox is now eliminated and eqs. ^ and (jU) can 
be solved analytically. 

For a supported membrane, Evans and Sackmann as- 
sumed that a thin lubricating layer of bulk fluid (thickness 
H and 3D viscosity ryt ) exists between the membrane and 
substrate [15] . In this case, they specifically identified the 
friction parameter as A — rjf/H provided that H is small 
enough. Here we consider that the parameter A represents 
all kinds of frictional interactions between the lipid head 
group and the surrounding fluid even for a free membrane, 
and hence has more general significance. Previously, the 
above Brinkman model was used to calculate velocity au- 
tocorrelation function of a disk [31 , diffusion coefficient of 
a polymer chain in membranes [32] , or dynamics of con- 
centration fluctuations in binary fluid membranes [33] . As 
pointed out in ref. [34], the translational symmetry along 
the membrane surface is broken for eq. ([3]) due to the 
friction term. We have implicitly assumed here that the 
velocity at infinity vanishes so that the friction term is 
proportional to v itself. We also note that eqs. §3§ and Q 
describe a purely 2D model so that the 3D hydrodynamic 
nature of the bulk fluid is not taken into account. 

As schematically presented in fig. [TJ we consider a cir- 
cular liquid domain of radius R and viscosity rf moving 
with a velocity — U = (— U, 0) in a thin membrane sheet 
of viscosity rj. As a result of the domain motion, a veloc- 
ity field is induced around it as well as inside the domain. 
Our purpose is to calculate the drag force experienced by 
the domain. To generalize our treatment, we allow for the 
situation where the momentum decay parameters are dif- 
ferent between the matrix (A) and the domain (A'). Here- 
after, quantities without prime correspond to those of the 
matrix, while quantities with prime refer to those of the 
domain. 

The fluid velocity can generally be expressed as a gra- 
dient of some potential tp and curl of some vector A = 
(0, 0, A) which has only a single component [3D]- Then the 
velocities in the matrix and the domain regions are ex- 
pressed as 

V = -V(y5 + VXA, (5) 



and 



-V<p' + V x A', 



(0) 
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Fig. 1. Schematic picture showing a section of an infinite mem- 
brane with viscosity r\ and momentum decay parameter A in 
which a circular liquid domain of radius R with viscosity 7/ and 
momentum decay parameter A' is moving. The liquid domain 
moves with a velocity —U in the a;-direction. 



respectively. Substituting eqs. fl5|) and ([5]) into eq. ((U), we 
obtain 



vV = o, vy = o, 



(7) 



which are the Laplace equations. With the use of eqs. ([3]) 
and (O, one can show that eq. ((3]) can be satisfied if the 
pressures are given by 



P = A<A p' = Ay, 



(8) 



while A and A' obey the following equations: 

(V 2 - k 2 )A = 0, (V 2 - k' 2 )A' = 0. (9) 

Here we have defined the (inverse) hydrodynamic screen- 
ing lengths for the matrix and the domain as ft = (A/ry) 1 / 2 
and k' = (A'/?/) 1 / 2 , respectively. We shall next seek for 
the solutions to eqs. (JT]) and ([9]) subject to the appropri- 
ate boundary conditions for the translational motion of 
the domain. 



3 Calculation of drag coefficient 
3.1 Velocity and stress tensor 

It is convenient to work in the cylindrical polar coordi- 
nates (r, 8) with the origin at the center of the circular 
domain. First we consider the matrix region where r > R. 
Under the condition that the velocity and pressure must 
approach zero at large distances, we write down the solu- 
tions to eqs. (jTj) and ([9]) as follows: 



<p = — cos ( 
r 



A = C 2 K 1 ( K r) sm9. 



(10) 



and tangential components of the velocity are given by 

COS0, (11) 

J, (12) 



V0 



where we have used the recursion relations among the 
modified Bessel functions (see Appendix A). Then the 
components of the stress tensor can be obtained as 

dv r 

a rr = —p+ 2?7— — 
or 



+r, ( l^ 1 ( Kr ) + ^ (1 ,„, 



a r g = T] 



I dv r dvg vg 
r 89 dr r 
~4Ci 



(13) 



2k 

+C 2 ( -^K^kt) + k^K^kt) + —K (Kr) 



sin 8. 
(14) 



Inside the liquid domain where r < R, the proper so- 
lutions to eqs. (|7|) and © subject to the condition that 
they do not diverge as r — > are 



<p' = C[r cos 8, A' = C' 2 h («V) sin f 



(15) 



Here, C[ and C 2 are unknown coefficients, and Ii(z) is 
the modified Bessel function of the first kind of order one. 
Since the corresponding radial and tangential components 
of the velocity are now 



C 2 



r 



cos 9, 



C' 



(16) 

sin 8, (17) 
the components of the stress tensor can be obtained as 



C[ - CWUn'r) + ^Ji(kV) 
r 



-v 



±I 1 ( K >r) + K ' 2 I 1 ( K 'r)-—I (K>r) 



cos 8, 
(18) 

sin^. 
(19) 



Here, C\ and C 2 are unknown coefficients which will be 
determined from the boundary conditions, and K\(z) is 
the modified Bessel function of the second kind of order 
one. Although the general solutions for cp and A can be 
expressed as series expansions in terms of r, we have kept 
only the least number of terms satisfying the requisite 
pressure and velocity conditions. From eq. (JSJ, the radial 



3.2 Boundary conditions 

Next we assume that the no-slip condition is satisfied at 
the domain boundary. This means that, at r = R, the 
radial component of the fluid velocities should be equal to 
the domain velocity — £/cos0, the tangential components 
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Fig. 2. The dimensionless drag coefficient £ as a function of 
dimensionless domain radius e = kR. The different curves are 
for E = 0.1 (solid), E = 1 (dashed), and £ = 10 (dotted). 



should be continuous, and so should the components of 
the stress tensor [53]. These conditions are written as 

v r = -Ucos6, (20) 

v' r = -Ucos6, (21) 

v e =v' e , (22) 

aye = a' r6 . (23) 

Since we have kept only the lowest order terms in r for 
the solutions to eqs. Q and ([S]) (see eqs. ([ID]) and (IT51 ). 
we are allowed to neglect the shape deformation of the 
circular domain while it moves. In other words, we are 
implicitly assuming that the line tension at the domain 
boundary exceeds a typical viscous force. A more quan- 
titative argument to justify this condition will be given 
later. 

Using the above four boundary conditions, we can ob- 
tain the four unknown coefficients C\ , C% , C[ and C 2 which 
are listed in Appendix B. In the following, we introduce 
the dimensionless parameters to measure the relative vis- 
cosities E = rj/rj' and the relative decay parameters L = 
A/ A'. We also use the notations e = kR and e' = k'R = 
ey E/L to measure the domain radius. Furthermore, the 
arguments of the modified Bessel functions will be omit- 
ted as K n = K n (e) and I n = I n {^) in order to keep the 
notation compact. 

3.3 Force exerted on a domain 

The force exerted on the domain is given by the integral 
of the stress tensor over the boundary: 

/•2tt 

F = R d9 {a„ cos B - o rd sin 9). (24) 
Jo 
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Fig. 3. The dimensionless drag coefficient ( as a function of 
the relative viscosity ratio E = rj/rj 1 when e = 1. The different 
curves are for L = 0.1 (solid), L = 1 (dashed) and L = 10 
(dotted). 

After some calculation, the drag coefficient obtained from 
C = F/U becomes 

{(e;E,L) = e 2 
Airr/ 4 

eifi[(4 + e' 2 )h - 2e% + 2E(e'I - 2JQ] 

+ K [(i + e' 2 )h - 2e'I ] + E(2K + eK x ){e'I Q - 2h) ' 

(25) 

An alternate form can be obtained by using the following 
recurrence relation 

W) = ^P-+W), (26) 

so that we have 

C(e;E,L) _ e 2 eK^e'h - 2J 2 + 2EI 2 ) 

4tt77 4 K {e'h - 27 2 ) + E(2K + eK x )h ' 

This is the main result of the paper. 

In fig. [2] we plot the dimensionless drag coefficient £ 
as a function of dimensionless domain size e = kR for 
E = 0.1, 1, 10 when L = 1. In all these cases, the drag co- 
efficient increases with the domain radius R as it should 
be. For fixed values of e, on the other hand, the drag co- 
efficient is smaller when the domain viscosity rj' becomes 
smaller (larger E). Fixing the domain size to e = 1, we 
have plotted in figs. [3] and Q] the drag coefficient ( as a 
function of E and L, respectively. In the former, we chose 
various values of L ranging from L = 0.1 to 10, while the 
values of E were changed from E — 0.1 to 10 in the lat- 
ter. Note that these graphs are both semi-log plots. From 
fig. [3] we see that the drag coefficient monotonically de- 
creases with increasing E (smaller domain viscosity) . This 
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Fig. 4. The dimensionless drag coefficient ( as a function of 
the relative momentum decay parameter ratio L = A/ A' when 
e — 1. The different curves are for E — 0.1 (solid), E = 1 
(dashed) and E = f (dotted). 
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Fig. 5. The limiting dimensionless drag coefficient Co (solid, 
eq. pop ) and Coo (dashed, eq. (131 [I 'l as a function of dimension- 
less domain radius e = kR when L = 1. 



can be attributed to the internal flows generated in the do- 
main because they are more efficient in transporting the 
fluid around it. As a result, a domain with finite viscosity 
feels a smaller drag force than a rigid disk. The asymptotic 
values of C for E — > and E — > oo converge to respective 
constants. In fig. [U the drag coefficient is a decreasing 
function of L. The large L values of C are dependent on 
E as discussed below. 



3.4 Limiting expressions 

Next we examine several asymptotic limits of eq. (|25l) or 
eq. ([27|) . Some useful formulae are listed in Appendix A. 

First we consider arbitrary values of E and L. In the 
limit of e — > 0, eq. ([2"5"1) gives a logarithmic behavior, 



C(e -> 0) 

47TJ7 



l + £ 



ln(2/e)-7 + £[ln(2/e)- 7 +(l/2)] 



. (28) 



where 7 = 0.5772 • • • is Euler's constant. The above in- 
dependent logarithmic behavior is a new outcome of our 
result. In the opposite e — > 00 limit, we have an algebraic 
dependence 



C(e -> 00) e 2 



4:777] 



4 ' 



(29) 



which does not depend on E. We also note that the above 
asymptotic expressions for small and large domain size 
limits do not depend on L. 

For arbitrary e, we next consider the rigid disk case 
(rf — !> 00) where the limit of E — > is taken. Then we 
obtain 



C(e;E^0) _ e 2 | eKi _ Co(e) 
Airrj 4 Kq Atttj 



(30) 



which coincides with the result by Evans and Sackmann 
[T§] . The opposite limit of E — > 00 corresponds to a 2D 
gas bubble (rf — > 0). In this case, we obtain a different 
expression: 



C(e;£^^) 

47T77 



2e^! 



_ Coo(e) 

4 2A' + eKi inrj ' 



(31) 



In fig. [5l we plot both Co and Coo as a function of e in log- 
log scales. In the limit of e — > 00, the difference between 
the two curves is the order of e which becomes negligibly 
small compared to the first terms in eqs. ([5D|) and (|3~lj) . 
Upon taking the e — > limit of eqs. (|3U)) and (|3"Tj) , we 
obtain 

1 (32) 



47T7/ 

Coo(e^O) 



ln(2/ C )- 7 ' 
1 



(33) 



Am] ln(2/e)-7+(l/2) 

These are nothing but the E — > and E — > 00 limits of 
eq. (|2^|) . We note that small e behavior of Co an d Coo do 
not coincide in fig. El The limiting expressions of eqs. (1301) 
and (j3"T|) for e — > 00 coincide with eq. (l29l) . 

In the limit of L — > 0, which is the case of high momen- 
tum dissipation in the domain region, we recover eq. (1301) 
again; 



C(e;L^0) 



K 



Atttj 4 
In the opposite L —¥ 00 limit, we obtain 

C(e;L 00) _ e 2 2e(l + E)K X 

4^ ~ T + 2(1 + £;)AT + eEKi 



which depends on E as observed in fig. [4] Equation ([35]) 
reduces to eqs. ([501) an d (ETJ) for £ -> and E -> 00, 
respectively. 



(34) 



(35) 
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4 Conclusion and discussion 

In summary, we have obtained the drag coefficient of a cir- 
cular liquid domain which has a finite viscosity using a 2D 
hydrodynamic equation with momentum decay or the 2D 
Brinkman equation. We showed that the drag coefficient 
decreases as the domain viscosity becomes smaller with 
respect to the matrix viscosity. This is because the flow in 
the domain helps to transport the fluid in the surrounding 
matrix more efficiently. Our result shows an independence 
in the small domain size limit (eq. ([28])). New limiting ex- 
pressions for the E — > oo (eq. (f3"Tj)) and L — > oo (eq. ([55]) ) 
are also obtained. 

Several points merit further discussion. We first discuss 
the realistic value of E — 77/77' in the case of lipid domains. 
Lipid domains are known to arise from a phase separation 
between a liquid-ordered (L Q ) phase rich in saturated lipid 
such as sphingomyelin and a liquid-disordered (Ld) phase 
rich in unsaturated lipid Applying the pulsed field 
gradient NMR spectroscopy, Oradd et al. measured lateral 
diffusion coefficients of a single lipid molecule both in the 
L a and Ld phases [7] . With the use of the SD logarithmic 
expression for the diffusion coefficient, the 2D viscosities 
of these phases are estimated as r) w 1.6 x 1CP 9 Ns/m and 
77d ~ 0.4 x 10~ 9 Ns/m at 293 K (the membrane thickness 
is chosen to be h « 3.8 x 10~ 9 m). When the L Q phase 
forms a domain in the matrix of the Ld phase, the typi- 
cal viscosity ratio would be E ss 0.2. In the opposite case 
where the Ld phase forms a domain, the ratio tends to 
be E « 4. Hence, it is necessary to take into account the 
finite viscosity ratio between the domain and the matrix 
rather than regarding the domain as a rigid disk. Inci- 
dentally, sphingomyelin constituting the L Q phases has a 
large head group which can protrude out into the 3D fluid, 
whereas the Ld phase are devoid of such structures [T2] . 
It is therefore reasonable to assume that the momentum 
decay parameter A experienced by the L Q and Ld domains 
are different. 

When imposing the boundary conditions as given from 
eq. (|2"0|) to (|23[) . we have assumed that the domain does 
not undergo any shape deformation. This implies that the 
line tension should be large enough to overcome the vis- 
cous force of deformation. In the case of lipid domains, 
the line tension was measured to be a ~ 10~ 12 N us- 
ing domain boundary flicker spectroscopy [35] . For a flow 
of the order of U ~ 10~ 6 m/s, a typical viscous force 
turns out be r/U ~ 10~ 15 N which is much smaller than 
a. Hence it is reasonable to assume that the line tension 
is large enough to maintain the circular shape of the do- 
mains unless the temperature is very close to the critical 
point. Different boundary conditions were used in order to 
consider the relaxation of deformed domains in a polymer 
monolayer [36| . 

The crossover from the logarithmic to algebraic be- 
haviors of the drag coefficient takes place when e ~ 1 or 
R ~ A typical domain size which corresponds to this 
condition is roughly estimated below. As stated in sect. [21 
the momentum decay parameter for a supported mem- 
brane was originally identified as A = rjf/H, where rjf is 
the 3D viscosity in the outer fluid and H is the thick- 



ness of a thin lubricating water which exists between the 
membrane and substrate [19] . In this situation, we have 
the hydrodynamic screening length k^ 1 — y/ j^H/r/f. Al- 
though we do not assume the presence of a solid substrate, 
the hydrodynamic screening length in the Brinkman ap- 
proach would typically be kT 1 ~ 10 -7 m when H = 10 -8 
m, 77f ~ 10 -3 Ns/m 2 and 77 ~ 10~ 9 Ns/m. Domain size in 
the order ofi?^10~ 7 mis accessible in the experiments. 

Within the present Brinkman model, the drag coeffi- 
cient has a R 2 dependence in the large size limit as in 
eq. (1291) . Such a behavior was indeed observed experimen- 
tally for Brownian motion of phase-separated domains on 
stacked-supported membranes [371 . However, this depen- 
dence differs from the asymptotic analysis of the SD model 
according to which the drag coefficient increases linearly 
with R (similar to a 3D sphere) [U]. For ternary vesicles, 
Cicuta et al. have measured the drag coefficient of a do- 
main which is proportional to R [5J, supporting the SD 
model. The large scale behavior of the Brinkman model 
and the SD model differs because the former is essentially 
a 2D model while the latter has a 3D character. It should 
be also noted that the hydrodynamic screening length in 
the SD theory is given by v^ 1 = i]/i]f which is different 
from kT x . In the present Brinkman approach, the hydro- 
dynamic screening length is the geometric mean of v^ 1 
and H. 

At this stage we mention that, based on the SD theory, 
an integral form of the drag coefficient for a 2D liquid 
domain with E = 1 was obtained by De Koker as [55] 



Cdk 

47T77 



dz 



z 2 (z + S) 



(36) 



where S = vR and J\ (z) is the Bessel function of the first 
kind. We have estimated the above integral numerically. In 
the limit of 5 —> 0, the above expression gives a logarithmic 
dependence on i?, in agreement with our result for small 
e as seen in eq. ([2"5]) . In the opposite limit of S —> 00, 
however, £dk is proportional to R due to the 3D nature 
of the SD theory. As describe above, this is different from 
eq. ([2"9"j) which gives the R 2 dependence. 

Several authors considered a more general scenario in 
which a disk or a domain moves in a membrane being lo- 
cated at a finite distance H from a fixed substrate [2(3 39140141142) . 
Among these, Stone and Ajdari numerically solved the 
coupled dual integral equations which are obtained as so- 
lutions for the Stokes equations [39^ • They showed that the 
Brinkman model gives a qualitatively good approximation 
over a wide range of membrane-substrate separations. One 
can explicitly show that the Brinkman and the SD mod- 
els are recovered in the limit of H — » and H — > 00, 
respectively. 

For comparison, we finally write down the analogous 
expression for the drag coefficient of a spherical drop hav- 
ing 3D viscosity 773 and radius R moving in a fluid of 3D 
viscosity 773 [23124129] 



C.3 



2L, 



27T77 3 i? E 3 + 1 ' 



(37) 
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where E 3 = r) 3 /r)' 3 . In the solid sphere limit (r)' 3 — > oo), 
we obtain (3 = Qnrj 3 R (the Stokes result). In the opposite 
bubble limit (r]' 3 — > 0), we have £3 = 4:irr] 3 R. The reduc- 
tion in the drag coefficient with the decrease in the drop 
viscosity rj' 3 is similar to that seen for a liquid domain in 
2D as discussed before. It is also important to point out 
that the effect of viscosity ratio E 3 exists for all the drop 
sizes R because there is no typical length scale for a 3D 
fluid. This is in sharp contrast to our model in which the 
hydrodynamic screening length k _1 plays a crucial role. 

We thank S. L. Keller and Y. Sakuma for useful discus- 
sions. This work was supported by KAKENHI (Grant-in- 
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Appendix A. Useful formulae 

|, the following recursion relations 
K' Q (z) = -K^z), 



According to ref. 
hold: 



K{(z) = -K (z) 
I&)=h(z), 

l[{z)=h[z) 



Hz) 



(A.I) 

(A.2) 
(A.3) 
(A.4) 



Here the primes indicate the derivative with respect to z. 
In the limit of z — > 0, we have 



K (z 

K x {z 

Hz 
h{z 
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" 7+ 4 


1 - 
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- 7 











z 2 + 



1 I 

~z ~ 2 



In 



z + 



1 



z 
1 
z z 
2 + 16 



9G 



(A.5) 
(A.6) 

(A.7) 
(A.8) 
(A.9) 



where 7 = 0.5772 ... is the Euler constant. In the limit of 
z — > 00, we have 



K n (z) 

In(z) 



'2irz 



1 

1 - 



An 2 - 1 

An 2 - 1 
~^z~ + 



(A.10) 
(ATI) 



Appendix B. Coefficients Ci,C 2 ,C[ and C' 2 

The coefficients C\, C 2 , C[ and C' 2 are given by 

Ci = -R 2 U - 2K l RU/{nK f) ) 

+ 2K 2 R 2 r 1 U(-2I 1 + 7 i?^)/[^o(M 1 + M 2 )], (B.l) 
C 2 = 4I Q R(r) - V ')k'U/[k(M 1 + M 2 )] 

- 2hU[A<n - ?/(4 + i? 2 K ' 2 )]/[ K (M! + M 2 )], (B.2) 

C[ = U + 2I 1 K 1 Rt 1 kU/(M 1 + M 2 ), (B.3) 

C 2 = 2K 1 R 2 r 1 nU/(M 1 + M 2 ), (B.4) 

with 

Mi = 2I Q K Q R(r] - t)')k' + hK^r^nn 1 , (B.5) 

M 2 = —2I\K\Rr\K - U x KqT] + hK r}'(4 + R 2 k' 2 ). 

(B.6) 

In the above, K n = K n (nR) ~ K n (e) and I n = I n (K'R) = 
/„(e') as denoted in the text. 
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